Manifestation of electron-electron interaction in the magnetoresistance of graphene 
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We investigate the magnetotransport in large area graphene Hall bars epitaxially grown on sili- 
con carbide. In the intermediate field regime between weak localization and Landau quantization 
the observed temperature-dependent parabolic magnetoresistivity (MR) is a manifestation of the 
electron-electron interaction (EEI) . We can consistently describe the data with a model for diffusive 
(magneto)transport that also includes magnetic-field dependent effects originating from ballistic 
time scales. We find an excellent agreement between the experimentally observed temperature de- 
pendence of MR and the theory of EEI in the diffusive regime. We can further assign a temperature- 
driven crossover to the reduction of the multiplet modes contributing to EEI from 7 to 3 due to 
intervalley scattering. In addition, we find a temperature independent ballistic contribution to the 
MR in classically strong magnetic fields. 

PACS numbers: 



Graphene is an outstanding electronic material with 
unique physical properties and with excellent prospects 
for technological applications [l|-[3|- The linear band 
structure and the chiral nature of charge carriers re- 
sult in a rich variety of phenomena which can be 
largely described within a model of non-interacting elec- 
trons. Evidence for electron-electron interaction (EEI) 
in graphene, however, is sparse. The most spectacular 
example is the occurrence of the fractional quantum Hall 
effect and the insulating phase in strong magnetic fields 
Tol 11 1 . Another line of research involves the investiga- 
tion of EEI-induced renormalization of the Dirac spec- 
trum in graphene [HI, [l]| . 

Furthermore EEI plays a prominent role in electron 
transport also in the absence of quantizing magnetic field, 
when the system is a Fermi liquid. In particular, EEI 
leads to a quantum contribution to zero-field resistivity, 
that in two dimensional electron gases (2DEG) at suffi- 
ciently low temperatures T displays a logarithmic T de- 
pendence . Logarithmic corrections may further stem 
from weak localization (WL), when observed in very low 
fields. At somewhat hi gher fields logarithmic p(T) data 
were interpreted as EEI |15l - [l7| or Kondo effect [l8( ■ This 
indicates that an evaluation of the temperature depen- 
dence p{T) alone is insufficient for identification of the 
underlying mechanisms. A way to unambiguously single 
out the EEI effect is to study the temperature-dependent 
magnetoresistance (MR) of a 2DEG in non-quantizing 
magnetic fields, between weak fields where WL is opera- 
tive and strong fields where Landau quantization is fully 
developed. The analysis of EEI effects in the magne- 
totransport, which have been widely disregarded, is the 
scope of this paper. 

For this investigation large sample dimensions arc 
mandatory, as in small flakes mcsoscopic conductance 



fluctuations would add to the MR and the finite sample 
size sets a low-energy cut-off to the EEI. This is the case 
for most experiments using exfoliated graphene. There- 
fore we use epitaxial graphene on insulating silicon car- 
bide (SiC) which is available with excellent homogeneity 
on the wafer scale Our special choice is quasi free- 



standing monolayer graphene (QFMLG) [20), |21| which 



represents a hole-conducting graphene material with a 
charge carrier density ofn~3- 10 12 cm -2 and a mobility 
fj, of approximately 2000 cm 2 /Vs. In contrast to other 
epitaxial graphene materials n and /i are more or less in- 
dependent of temperature due to a reduced coupling to 
the substrate 2l| . The chiral nature of charge carriers in 
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QFMLG has been proven in previous experiments 

For the present study we prepared Hall bars by stan- 
dard electron beam lithography Q and measured the 
magnetoresistivity p(B) in four-point geometry for var- 
ious temperatures T using low-frequency lock-in tech- 
nique. For all p(B,T) curves the magnetic field B was 
swept from -8 to 8 T twice to confirm reproducibility. 
All measurements were performed under cryogenic vac- 
uum in order to reduce the effect of adsorbatcs. In total 
we investigated four samples of different size with very 
similar findings as concerned the EEI effects [24[ . In this 
paper we focus on one particular sample SI (330 ^m long 
and 50 /im wide channel) which shows all observed phe- 
nomena in one consistent data set. Figure [TJa) shows 
the magnetoresistivity raw data. Obviously the Lan- 
dau quantization regime is not yet reached (however, it 
has been observed in higher fields (22|). This indicates 
that the damping of Shubnikov-de Haas oscillations in 
the sample is dominated by smooth disorder, so that 
the transport relaxation time is longer than the quan- 
tum scattering time. 

At weak fields the WL correction is clearly visible as a 
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zero-field peak. We have characterized it carefully and 
can accurately describe the anomaly within the stan- 
dard graphene model [1, [25j . The evaluation shown in 
Fig. &b) yields the intravalley scattering time responsi- 
ble for the breaking of the effective time-reversal sym- 
metry within a single valley = 0.075 ps and the in- 
tervalley scattering time t- iv = 1.1 ps, together with the 
phase relaxation time t v , which is strongly temperature 
dependent [cf. Fig. E^c)] and decays with increasing T 
roughly as T -1 as expected for electron-electron scatter- 
ing in the diffusive regime [1 , 14 1 . The WL data are qual- 
itatively consistent with findings in exfoliated graphene 
where the typical time scales are r lw = 0.8 ... 2.5 ps and 
r* = 0.05 ... 0.3 ps [1, El. The deviation from the T" 1 
dependence at low T presumably stems from the presence 
of magnetic impurities fl7| . 

At stronger fields we observe that the curves have 
roughly parabolic shape and vary with temperature, 
showing an evolution from negative to positive magne- 
toresistivity as the temperature rises. Thus, the observed 
InT dependence (see the analysis below) indicates that 
we are dealing with the quantum contribution due to 
EEI as classical contributions to the magnetoresistancc 
are temperature independent at low T . Another striking 
feature of Fig. QJa) is a crossing point of low-T curves 
that moves to a higher value of magnetic field at higher 
temperatures. 

For further analysis the data are symmetrized with re- 
spect to B and relative differences (p — po)/Po ~ — & a 
are considered. We choose po = p(0.5T) to eliminate 
the effect of the WL correction at \B\ < 0.5 T. The re- 
sulting curves (p — po)/po(B) are plotted in Fig. QJb). 
Again the transition from negative to positive magne- 
toresistance with rising temperature is visible. A closer 
look at intermediate temperatures reveals that the curva- 
ture differs in its sign even within the same curve between 
weaker (\B\ < 4.5 T) and stronger fields. Indeed for all 
curves two separate parabolas (w c T tr ) 2 ■ M with dif- 
ferent curvatures M\ and M Q have to be taken into ac- 
count for the inner and the outer field region respectively, 
where w c = j^=B is the cyclotron frequency. As a par- 
ticular example Fig. [TJc) shows this decomposition for 
the lowest (T = 1.8 K) and an intermediate (T = 20 K) 
temperature. The transition field of ~ 4.5 T is constant 
for all temperatures. 

We analyze the data with the help of a model that 
successfully described the magnetoresistance in many 
2DEGs in non-quantizing fields (see Refs 27 1 and [28| for 
review) and is equally suited for graphene at large dimen- 
sionless conductance. The model takes EEI into account 
and covers the crossover from the diffusive to the ballistic 
regime. The diffusive regime where k^TT^/h <C 1 is ob- 
served for the lowest-temperature curves. The transport 
time Ttr ~ 0.045 ps is derived from the Drude resistivity. 
The highcr-T curves correspond to the crossover regime, 




FIG. 1: Evolution of magnetoresistivity with temperature, 
(a) The curvature changes sign with increasing temperature. 
Further, for temperatures T < 10 K the curves cross at B « 
5 T; for higher temperatures the crossing point shifts towards 
a larger value B ~ 9 T. (b) Relative magnetoresistivity (p — 

Po)/ 'Po- ( c ) Two parabolas (aj c Tt r are fitted to the 

data in (b) (shown representative for 1.8 K and 20 K): one in 
the diffusive regime for weak fields 1 T < \B\ < 4.5 T (dashed) 
and one in the diffusive regime with ballistic contributions 
for \B\ > 5T (dotted). Two sketches illustrate the different 
scattering mechanisms in both regimes. 



where the importance of ballistic effects [27, in- 
creases. Still, for all studied temperatures and samples 
the condition k^TT^/h < 1 is fulfilled which means that 
the ballistic regime is not reached. 

In the diffusive regime, EEI induces a logarithmic- 
in-T, B-indepcndent correction <5<7 ee to the longitudinal 
conductivity, while not affecting the Hall conductivity 



3 



fl4| . As a result, EEI gives rise to a magnetoresistance 
Sp/po oc (p,B) 2 S<T cc /<To (here <tq is the zero-f? Drude con- 
ductivity) which is parabolic in B [32], |33j : 



p = po + [(w c r tr ) 2 - 1] 



A -In 



(1) 



This correction includes two contributions: the Coulomb 
term and the Hartree term. The coefficient A is the sum 
of these two terms 



A = A c +A H = l + c 



1 



ln(l 



(2) 



with the Fermi liquid constant Fq . It is a measure for the 
interaction strength and in the case of weak interaction 
it can be calculated within Thomas-Fermi approxima- 
tion by integrating the interaction matrix element over 
all scattering angles [l6| : 



F t 



0, thco 



-a 



'e/2 



2tt (sin 6>/2 + 2a) 



AO, 



(3) 



with the interaction constant a = e*/(4n£oHvF) ~ 0.41 
and an effective charge e» = ey2/(e + 1) w 0.43 e tak- 
ing into account the screening of Coulomb interaction 
by the SiC substrate (esic = 9-66). For the given case 
equations © and © yield A(F° thoo ) w 1 - 0.046 • c 
and Fq theo w —0.09 respectively. The constant c is the 
number of multiplet channels participating to EEI. 

Here the peculiarity of graphene enters at the follow- 
ing two points. First, the cos 2 0/2 reflects the influence 
of chirality on scattering of charge carriers in graphene. 
Second, the channel statistics of graphene is different 
from a standard 2DEG and depends on the interplay of 
symmetry-breaking time scales. In our samples the hier- 
archy is Tt r < t* < Tj v as extracted from the analysis of 
the WL peak. This hierarchy is also found in exfoliated 
graphene 0. In perfect graphene the spin and valley 
degeneracy gives rise to 15 multiplet channels. For tem- 
peratures T < hjk^T-,, the breaking of the time-reversal 
symmetry within a single valley (due to, for example, the 
trigonal warping or ripples) suppresses 8 channels out of 
15, yielding c — 7 [16| . The intervalley scattering fur- 
ther reduces the number of multiplet channels to c = 3 
(spin-triplet). This standard 2DEG case is expected for 
sufficiently low temperatures T < T- lv := h/k-QT^ ~ 7K. 

When analyzing the prefactors A of formula ([1]), one 
can expect to see the predicted crossover between c = 3 
and c = 7 in the temperature dependence of the curva- 
ture of the magnetoresistivity. Figure [5] plots the cur- 
vatures Mi for the inner (circles) and M Q for the outer 
(squares) field region as a function of temperature (more 
explicitly k^TT^/h). In this plot two regions of different 
slope can be distinguished, with a crossover at T lv « 7K, 
that perfectly fits the temperature scale given by r; v (ex- 
tracted from WL data). From Eq. (TTJ it is obvious that 
the slope of the curvatures M in this plot over k^Tr^/h 
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FIG. 2: (a) The curvature M of the parabolic MR extracted 
from the fits to the curves shown in Fig. \V[b) at weak (circles) 
and strong (squares) magnetic field. A kink is observed at 
T — ft/rivfca ~ 7K as the number of contributing multiplet 
channels changes from c = 3 to c = 7. From the slopes A 
the Fermi-liquid constant Fq is calculated for the different 
regions, (b) Weak localization correction to the conductivity 
at weak fields. Straight lines are fits of the graphene-specific 
model [25| that allows to determine the scattering times t*, 
Tiv and the phase-coherence time t v . (c) As expected for 
electron-electron scattering in the diffusive regime t v decays 
as T~\ 



(on logarithmic scale) is equal to the coefficient A. Hence 
the change in slope is caused by a different number of 
contributing multiplet channels in the two regimes. Us- 
ing c = 3 as expected for T < T lv a Fermi liquid constant 
of Fqj = -0.05 and F a o = -0.06 can be extracted for 
the inner and outer field region respectively. For T > T; v 
F CT ; = -0.14 and F CT o = -0.16 are found with c = 7. 
The obtained values are reasonably close to the value of 



the theoretically expected FQ thco = —0.09. The higher- 
values at higher temperatures can qualitatively be at- 
tributed to the crossover towards ballistics (see discussion 
below). 

It is worthwhile to compare the extracted values of Fq 
to those measured in a 2D electron gas in conventional 
semiconductor structures. The main factor determining 
the Fermi liquid constant is the widely used parameter 
r s controlling the ratio of characteristic potential and ki- 
netic energies. It depends on the electron charge, the 
dielectric constant, and the Fermi velocity (in our case 
r s = \[2a). The central difference to semiconductors is 
that there the Fermi velocity depends strongly on the 
electron concentration n, so that r s varies substantially 
(typically in the range from 2 to 20, leading to a change 
of Fq between 0.15 and 0.7; see, e.g., Ref. [34|). In con- 
trast, in graphene the velocity is roughly n- independent, 
and r s remains relatively small. In addition, for given 
r s the Fermi-liquid constant is somewhat suppressed in 
graphene due to chiral character of carriers, see Eq. ([3]). 
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A hallmark of the diffusive regime is that the curva- 
ture of the MR correlates with the sign and amplitude 
of the correction to the zero--B resistivity. In particular 
this leads to the well defined crossing point (CP) in Fig. 
□[a). At this point at B C p = 5.2 T [Fig. EJa)] where In- 
dependent and B-dependent corrections exactly cancel, 



3J. The 



the resistivity has the classical Drude value [3 
value of this CP agrees well with Bcp = 4.8 T expected 
from cj c T tr = 1. The overall consistency of our find- 
ings (together with the emergence of the crossing point) 
confirms that the applied diffusive model describes the 
physics properly and hence is a reliable way to separate 
the effect of EEI from other corrections. 

Now we turn to the observations which are beyond the 
purely diffusive model described by Eq. (JTJ) . (I) First, we 
find different curvatures of the MR for weak and strong 
magnetic fields [cf. Fig. QJc)]. (II) Second, the higher- T 
MR curves have a tendency to form an apparent second 
crossing point [cf. Fig.|3jb)]. These two observations sug- 
gest that the diffusive mechanism, which governs the log- 
arithmic temperature dependence of the curvature of the 
MR in the whole range of magnetic fields and tempera- 
tures explored in this work, is supplemented by ballistic 
effects that become more pronounced at higher T and 
magnetic field. 

Observation (I) implies that the MR contains a T- 
independent contribution which stems from ballistic 
scales 
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3l| and is sensitive to the relation between 
uj c and Tt r . As the magnetic field increases, for B > 
(which corresponds to w c T tr > 1) electrons can accom- 
plish closed cyclotron orbits which leads to an enhanced 
probability of scattering angles 9 « [sketch in Fig.QTc)]. 
This correlation causes a change of the strength of EEI 
[see Eq. ([3])] in the ballistic contribution and hence the 
overall curvature of the MR. In addition, there might be a 
temperature-independent contribution due to a classical 
mechanism of magnetoresistance (due to non-Markovian 
processes, for review see Ref. [H]) operative in classically 
strong magnetic fields B > that leads to a vertical 
shift between the values of Mj and M a . 

Observation (II) indicates that the T dependence of 
the zero-B resistivity on does not fully determine the 



curvature of the MR |27J, |31|, |36|, |37| . This signals the 



onset of a crossover to the ballistic regime at T > 20 K. 
Indeed, at k^TT^jfi < 1 the T dependence of po stems 
from the combination of a logarithmic diffusive correc- 
tion and a T-dependent correction Sp(T) to the mobility 
(Ballistic return processes may give rise to ballistic In- 
dependent corrections to the mobility (38j). At lowest 
temperatures, the logarithmic diffusive correction dom- 
inates. With increasing T, the diffusive term decreases 
while Sfi(T) grows. Importantly, the processes giving rise 
to a T-dependent correction to the mobility do not lead 
to the MR at finite B, in contrast to the logarithmic dif- 
fusive term. Therefore, 5p(T) only shifts all resistivity 
curves p(B) in the vertical direction (in the crossover 
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FIG. 3: Zoom on the crossing points, (a) The crossing point 
at Scp = 5.2 T is already visible in Fig. \T\ (b) For higher 
temperatures, an apparent crossing point at Bcp2 — 9. IT 
emerges by extrapolating (dashed lines) the parabolic shape 
of the curves onto B > 8 T. 



to the ballistic regime the diffusive and ballistic con- 
tributions to this shift become comparable), while the 
T dependence of the MR curvature is still given solely 
by the logarithmic term. This explains the shift of the 
CP towards an apparent high-T CP at Bqp2 — 9.1 T 
that emerges when the MR parabolas are extrapolated 
to stronger fields [Fig. [3](b)] . Its observation thus proves 
that ballistic corrections to po gain importance at higher 
temperatures. 

Theoretical results for MR in the ballistic (and 
crossover) regimes depend on the microscopic model of 



disorder [27], [28[. Therefore, a systematic study of MR 
in the crossover and ballistic regimes may be helpful for 
clarifying the character of disorder in the system (which is 
particularly important for electron transport in graphene 
3^|): we relegate this to a future work. In this paper, we 
restrict ourselves to applying the diffusive model supple- 
mented at strong B by a T-indcpcndcnt ballistic contri- 
bution, which is a reasonably good approximation in the 
considered temperature range. 

In conclusion, we have studied the magnetotransport 
in graphene epitaxially grown on silicon carbide. We ob- 
serve a temperature-dependent parabolic magnetoresis- 
tivity which is a manifestation of the EEI. The tempera- 
ture dependence of the MR curvature is logarithmic, as 
expected for the diffusive regime. We find a pronounced 
kink in this dependence, that is related to the reduction 
of the number of multiplct modes contributing to EEI 
from 7 to 3 due to intervalley scattering. 

We show that our model of diffusive transport with bal- 
listic contributions can describe the underlying physics 
properly. The analysis of the temperature dependence of 
the magnetoresistance allows for an accurate description 
of the EEI interaction correction in the diffusive regime. 
At magnetic fields higher than the observed crossing 
point (uj c Ttr > 1) the magnetoresistance shows an addi- 
tional (compared to weak field) T-independent parabolic 
contribution due to ballistic effects. Further quantitative 
analysis of the crossover and ballistic regimes should be 
helpful for shedding light on the nature of disorder in 
graphene structures. 
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Our work demonstrates that the analysis of MR in the 
intermediate range of B is a convenient tool for explor- 
ing many-body effects in graphene. In particular, this 
method allows to precisely distinguish the effect of EEI 
from additional corrections to the temperature depen- 
dent resistivity that also have logarithmic temperature 
dependence. Furthermore, our approach allowed us to 
determine the Coulomb interaction strength in our sam- 
ples with a rather high accuracy. The value of the in- 
teraction constant obtained from the fit to the experi- 
mental data is F£ —0.06 at low T and —0.15 at high 
T (where the system approaches the ballistic regime), 
which is close to the theoretical estimate for graphene on 
silicon carbide, Fq « —0.09. These findings should also 
directly apply to other graphene samples (e.g., to exfoli- 
ated graphene) provided the system sizes are sufficiently 
large. 
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Supplemental Material for 'Manifestation of 
electron-electron interaction in the 
magnetoresistance of graphene' 

This supplemental material is dedicated to elucidate 
the generality and reproducibility of the findings reported 
in the main text. We have observed the parabolic shape 
of the magnetoresistance (MR) in between the weak- 
localization (WL) regime and the Landau quantization 
regime and its logarithmic temperature dependence for 
many samples. We have investigated it thoroughly for 
four samples. We decided to show data of sample SI as 
one consistent data set in the main text. The data of 
the remaining three samples (S2, S3 and S4) are subject 
to this supplemental material. The quasi free-standing 
monolayer graphene (QFMLG) sample S2 had smaller 
size than SI which made the analysis of the low tempera- 
ture behavior difficult and increased the amplitude of the 
temperature-independent quasiclassical MR. The mono- 
layer samples S3 and S4 (both not QFMLG) showed a 
large phonon contribution to the MR in view of a strong 
coupling to the substrate. Therefore, sample SI was 
best suited for our analysis of electron-electron interac- 
tion (EEI) effects. Nevertheless, samples S2, S3, and S4 
show the EEI effects in the temperature dependence of 
the MR quantitatively consistent with those in sample 
SI (details are presented below). 



Similar sample type 

Figure [4] displays results of QFMLG sample S2 that 
is similar to sample SI presented in the main text. The 
main difference is that S2 is smaller (10 ^m wide, 70 (mi 
long) which made it difficult to analyze the low temper- 
ature behavior and increased the amplitude of the T- 
independent quasiclassical MR. In particular at lowest 
temperatures universal conductance fluctuations become 
apparent. Along with the somewhat lower mobility, the 
crossing point has moved to higher fields, and a field 
driven cross-over can not be detected. The evaluation of 
the curvature of the MR, however, shows the same cross- 
over from c = 3 to c = 7 with rising temperature. The 
analysis yields similar values of the Fermi-liquid constant 
Ffi , at least for high temperatures. At low temperatures 
the data quality and density is unfortunately not as good 
as for sample SI and therefore does not allow to extract 
Fq for the c = 3 region. 



Different epitaxial graphene material 

Beyond QFMLG samples SI and S2 presented in the 
main text and above, we investigated two samples S3 and 
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FIG. 4: (a) The magnetoresistivity raw data exhibits a 
parabolic shape and an apparent crossing point at ~ 9T. (b) 
The relative magnetoresistivity shows the evolution of the cur- 
vatures with temperature, (c) Evaluation of the curvatures M 
yields the Fermi-liquid constant. Intervalley scattering time 
Ti v = 1.5 ps and intravalley scattering time r* = 0.051 ps were 
extracted from WL data. 

S4 using monolayer graphene (MLG) [l9| . Microscopi- 
cally the materials differ in an interface layer underneath 
the graphene [2(| Hlf . Phenomenologically, they differ in 
the sign of charge carriers (holes for QFMLG and elec- 
trons for MLG), and in different coupling to substrate 
phonons 0. As a consequence of the latter, MLG has 
a substantial dependence of the resistivity p on temper- 
ature T. It can readily be seen [Fig. EJa)] that a strong 
temperature dependent offset of the MR curves result. 
Therefore the occurrence of a crossing point is obscured. 
Figurcs[5]Jb) and (d) show the further analysis carried out 
in exactly the same manner as presented in the main text. 
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FIG. 5: The magnetoresistivity raw data of sample S3 (a) and S4 shows a pronounced vertical shift due to the coupling to 
substrate phonons in addition to the parabolic MR. At low temperatures and high fields Shubnikov-de Haas oscillations are 
visible. The relative magnetoresistivity (b) displays the evolution of the MR curvature with temperature. As MLG samples S3 
and S4 are processed on the same chip their electronic properties are virtually identical as illustrated in (c) where MR curves 
for 10 K, 1.5 K, 100 K and 200 K (from top to bottom) are compared. Thus, also the evaluation of the MR curvatures M of 
samples S3 and S4 are nearly equal (d). Note that the extracted values are very similar to the findings on sample SI. The 
scattering times t;„ = 0.53 ps and r* = 0.025 ps were extracted from WL data. 



The enhanced mobility has two main consequences: (i) 
quantum oscillations occur at substantially lower fields 
as in SI. (ii) the obscured cross-over point also moves 
to lower fields and therefore the field range of the inner 
parabola (between WL regime and the obscured crossing- 
point) is compressed. Hence it is difficult to evaluate the 
inner curvature Mi and we restrict ourselves to M a . The 
evaluation of the Fermi-liquid constant results in excel- 
lent agreement with the values extracted for SI. This dis- 
plays the robustness of the findings taking into account 
that a different graphene material, a different charge car- 
rier density and a different charge carrier type is present. 
The downturn in Fig.[5][d) happens close to ksTTtr/h = 1 
where the diffusive picture becomes inadequate. In addi- 
tion, the phonon effects at high temperatures may affect 



the T dependence of the MR curvature. 



Reproducibility 

One advantage of epitaxial graphene is the ability to 
very well reproduce data. As an example we show the 
MR of two samples S3 and S4 with nominally identi- 
cal geometry (50 /im wide, 330 /mi long) being processed 
on the same chip. Differences can hardly be recognized 
[Fig. [5^c) and (d)]. The above examples, however, show 
that when different fabrication parameters or different 
processing are chosen, the resulting samples may differ, 
most evidently in mobility and/or charge density. 



